A Self Healing Model Based on Polymer-Mediated Chromophore Correlations 
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Here we present a model of self healing in which correlations between chromophores, as mediated 
by the polymer, are key to the recovery process. Our model determines the size distribution of the 
correlation volume using a grand canonical ensemble through a free energy advantage parameter. 
Choosing a healing rate that is proportional to the number of undamaged molecules in a correlated 
region, and a decay rate proportional to the intensity normalized to the correlation volume, the 
ensemble average is shown to correctly predict decay and recovery of the population of disperse 
orange 11-DOll (l-amino-2-methylanthraquinone) molecules doped in PMMA polymer as a function 
of time and concentration as measured with amplified spontaneous emission and linear absorption 
spectroscopy using only three parameters that apply to the full set of data. Our model also predicts 
the temperature dependence of the process. One set of parameters should be characteristic of a 
particular polymer and dopant chromophore combination. Thus, use of the model in determining 
these parameters for various materials systems should provide the data needed to test fundamental 
models of the underlying mechanism responsible for self healing. 

PACS numbers: 



o 

> 
00 

p 
o 



X 



I. INTRODUCTION 

Photodegradation is an inherent problem of materials 
that are used for high-light-intensity applications such 
as organic lasing media [ll |2| or materials for use in 
nonlinear-optical devices. [3, 3 Peng first observed the 
recovery of fluorescence after photodegradation of a dye- 
doped polymer optical fiber. 5] Howell observed full re- 
covery after photodegradation of Amplified Spontaneous 
Emission (ASE) from Disperse Orange 11 (DOll) dye 
doped in PMMA polymer.|i, Q The DOll molecule was 
shown not to recover in liquid solution. [8] An intrigu- 
ing demonstration that laser cycling could make materi- 
als more efficient and robust to future degradation^, 
motivated further research. 

The mechanism responsible for the decay and recovery 
process is still not well understood. [Tlj The purpose of 
this contribution is to develop a model based on physi- 
cally realistic assumptions that take into account all ob- 
servations. The model that we propose correctly pre- 
dicts all observations, makes new predictions that can be 
tested, and is expressed in terms of only three parameters 
that may be calculable from first principles. 



II. BACKGROUND 

ASE is the most sensitive probe of photodegrada- 
tion and self healing that has been extensively used 
to characterize PMMA polymer doped with DOll 
chromophores. [llj While not as well studied, decay and 
recovery of the AF455 chromophore has been observed 
using two-photon absorption (TPA) spectroscopy, d, 
Since ASE and two photon absorption are nonlinear- 
optical processes, this makes them sensitive probes but 



can also lead to larger experimental uncertainties. Addi- 
tionally, identically-prepared polymers tend to vary from 
sample to sample, so, it may be difficult to reproduce 
data runs with a high degree of precision. 

Keeping these difficulties in mind, it is nevertheless 
possible to determine common features over the full set 
of observations that hint at the mechanisms responsible. 
In particular, the ASE studies suggest the following, 

1. Dye molecules that irreversibly photodegrade in a 
liquid are found to self heal after photodegradation 
in a polymer. 0, Q 

2. Full recovery of the ASE signal after self healing is 
observed only for dye concentrations near the satu- 
ration limit in the polymer. The degree of recovery 
decreases at lower concentrations. 

3. For a particular concentration, the decay rate de- 
pends more or less linearly on the pump intensity 
but the recovery rate is a constant. [ll| 

4. The two-photon absorption cross-section of AF455 
decays at a rate in proportion to the intensity, but 
recovers at a single recovery rateprovided that the 
sample is not severely damaged 

5. Linear dichroism is constant during decay and re- 
covery, suggesting that molecular reorientation is 
not responsible. |11| 

6. Decay of ASE signal is accompanied by a change in 
the absorption spectrum with an isosbestic point, 
suggesting that the decay product is a different 
species, or a perturbed version of the chromophore, 
and rules out the possibility that dye diffusion is 
responsible[ll[ unlike in some other svstemsfl^. 
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7. Evolution of the spatial population profile of a burn 
mark during recovery does not fit the diffusion 
model - additional evidence that diffusion is not 
the cause of recovery, [l^ 

8. At high-enough intensities that result in visible 
damage (burn marks and laser ablation), the sam- 
ple appears to decay irreversibly, but there is some 
evidence that the recovery times in these cases are 
much longer than for the regime in which the sam- 
ple degrades by a small amount, [l^. [Tsj 

9. After cycles of photodegradation and recovery, the 
ASE intensity increases and the decay time con- 
stant increases relative to a non-cycled sample. 

10. As a function of temperature, the ASE intensity de- 
creases, but, the change in linear absorption spec- 
trum peaks at a different wavelength than when the 
material is photodamaged, suggesting that a new 
species is populated at higher temperature (as we 
will show below). Furthermore, the species found 
at elevated temperature recovers more quickly than 
a photodamaged molecule. Thus, the thermally- 
populated species - which we call a bystander state, 
is not a damaged molecule. 

11. The decay rate is observed to decrease with increas- 
ing concentration, as we discus later (and shown in 
Figure [TT|l : but, the recovery rate increases with 
concentration. 

12. The change in the linear absorbance spectrum with 
concentration peaks at the same wavelength as does 
the change in the spectrum as the material photo- 
damages. No such change in the normalized linear 
absorbance spectrum is observed as a function of 
concentration in MMA liquid as shown in Figure 
HJ This suggests that the damaged species and dye 
aggregates may be related. 

The key observations that drive our model are that 
(1) full recovery requires high chromophore concentra- 
tion (Item [2] above), (2) the decay rate decrease with 
concentration (Item [TT] above) and (3) optical absorp- 
tion spectrum evolves as a function of time in the same 
way that it changes with concentration (Item [12] above). 
These observations, along with the fact that the poly- 
mer plays a critical role, suggests that the interactions 
between chromophores as mediated by the polymer are 
responsible for self healing. As such, we call our new 
model the correlated chromophore model, which is a gen- 
eralization of the simpler embedded chromophore model 
used in the past.[ll| The role of temperature will be 
taken into account using the grand canonical partition 
function, which includes correlation scales that depend 
on interactions between chromophores. 
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FIG. 1: (a) Change in absorbance as a function of dye con- 
centration with respect to the lowest DO 11 dye concentration 
of DOll dissolved in the liquid monomer MMA. (b) Same as 
(a) but for DOll doped in PMMA polymer, (c) Change in 
the absorbance due to photodegradation in DOll/PMMA. 



III. THE MODEL 

In the following sections, we begin by describing the 
embedded chromophore model, generalize it to include 
correlations between chromophores, apply the grand 
canonical ensemble to get the temperature-dependent 
distribution of the size of the correlated regions, and com- 
bine these into the full theory of the decay and healing 
process. 



A. Embedded Chromophore Model 

Consider a system of N non-interacting molecules, n 
of which are undamaged. Assuming that there is only 
one damaged species, the damaged population is N — n. 
If the recovery rate is /3 and the decay rate proportional 
to the intensity, /, and given by al, then the evolution 
of the population is given by. 



dn 
Iti 



-ain + I3{N -n), 



(1) 



3 



which has a solution for / 7^ of 

where the sample is assumed to be originally pristine at 
t = 0. 

The population of undamaged molecules, starting 
when the pump is turned off at time t^, is given by 

^ = l-(l-^)e-^(-o), (3) 

where n(to) is the initial undamaged population at time 
to- This model is consistent with the time evolution of the 
ASE intensity in DO 11 [11] and the TPA crossection of 
AF455 9] chromophores at fixed concentration and tem- 
perature. Even though the population model is by ne- 
cessity dependent on more than two populations (i.e. ex- 
cited electronic and vibronic states of each species must 
play a role), if optical excitations and de-excitations are 
fast compared with the photodegradation and recovery 
process, the dynamics of these other states can be ig- 
nored. Note that the population of the chromophore and 
the decay product can be differentiated by their unique 
optical properties, such as ASE efficiency and TPA cross 
section. 



B. Correlated Chromophore Model 

The noninteracting model given by Equation [1] can be 
phenomenologically generalized by allowing the rates a 
and /? to depend on intensity - for example, by expressing 
them as a series in the intensity if higher order contribu- 
tions are small, a is independent of intensity if the optical 
damage process is dominated by one-photon absorption. 
If two-photon absorption is important, than a correction 
term linear in the intensity must be added to the theory. 
Similarly, it is possible that the recovery rate is acceler- 
ated or suppressed in the presence of light, in which case 
an intensity-dependent correction term needs to be added 
to /3. The data for DOll/PMMA and AF455/PMMA at 
fixed concentration and temperature agree with the the- 
ory for constant a and /?, so the higher-order correction 
terms are ignored in this paper, but can be easily gener- 
alized if needed. 

The observation that decay and recovery kinetics de- 
pend on the concentration can be taken into account by 
making a and /? a function of the concentration of the 
chromophores and their decay products. In particular. 
Figure [T] shows that the absorption spectrum changes as 
a function of DOll concentration in PMMA polymer, an 
indication that the DOll molecules are interacting with 
each other. The fact that higher concentration samples 
show accelerated recovery fll] and decreased decay rates 
suggests that interactions between molecules may be re- 
sponsible for self healing. The nature of this interaction 
is not important in building a phenomenological model 




FIG. 2: Domains of dye in polymer. 



that depends on concentration and characteristic ener- 
gies. Indeed, the phenomenological parameters that we 
will define in our model would be calculable from first 
principles as tests of the nature of these interactions. 

In the generalized model, we redefine N to be the num- 
ber of molecules that are associated with each other, 
be it through physical aggregation of DOll molecules 
into dimers or microcrystalites, a correlated region of 
dyes that interact through the polymer as mediated by 
phonons, locally oriented domains, or altogether new 
physics. For the purposes of this paper, we will generi- 
cally call such a correlated region a domain. Each domain 
will have its own characteristic decay and recovery rate 
depending solely on the domain size; and, the domain 
size distribution will be determined by the grand canon- 
ical ensemble. The observed bulk behavior will then be 
given by the ensemble average. Figure [5] shows a car- 
toon representation of a collection of domains with one 
of the domains specified with N molecules, n undamaged 
molecules and N — n damaged species. 

First, we focus on a single domain of fixed size N. We 
propose that the recovery rate of a damaged species will 
be accelerated in the presence of undamaged molecules 
in proportion to the number of undamaged molecules. 
Thus, in a bigger domain, the recovery rate will be larger 
than in a small domain. On the other hand, a small 
domain with mostly undamaged molecules will recover 
at a faster rate than larger domains that are populated 
with a preponderance of degraded molecules. 

There are several mechanisms that could lead to such 
behavior. For example, if the undamaged molecules are 
strongly interacting with each other, forming a damaged 
species will cost energy. The more neighbors, the higher 
the energy cost, leading to a slower decay rate and faster 
recovery rate. The nature of these interactions is not 
yet understood, but there are many potential candidates, 
from electric forces that are typically responsible for ag- 
gregation into dimers and microcrystalites to spin statis- 
tics that cause Bose-Einstein condensation or perhaps 
new unknown mechanisms. Based on experimental ob- 
servations, we find that the decay rate a/, depends in- 
versely on the number of molecules in a domain, or that 
the decay rate is given by al/N. These two generaliza- 
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tion of Equation [T] leads to 

dn ^ s al , , , 

— = l3n{N -n)- —n. (4) 

Integrating Equation |4l yields 

(N - alipN) no 

n 

710 + (iV - no - al/PN) exp [- {N P - al/N) t] ' 

(5) 

where rig is the initial undamaged population at i = 0. 
When the pump is turned off (/ = 0), Equation [S] ap- 
proaches n = at infinite time. With the pump on, 
at the infinite time limit, the population n approaches 
N — al / (3N provided that the intensity is below the crit- 
ical intensity I < Ic = N'^fi/a and vanishes at large time 
when the intensity is above this critical value. Past work 
showed that for I < Ic, the ASE intensity for DO 11 in 
PMMA polymer did indeed asymptotically approach a 
nonzero value at long times. [ll| 

With / = 0, the population recovers according to 



— = j3n{N - n). 
at 



(6) 



With n{to) as the undamaged population after the decay, 

Nn{to) 



n{to) +[N- n{to)] e(-^'3t) ■ 



N 



N 



n{to) 



4-NI3t) 



(7) 



(8) 



C. Bystander state 

A plot of the optical absorbance as a function of tem- 
perature, as shown in Figure shows the main absorp- 
tion peak decreasing as a new peak grows in proportion 
to the main peak's decrease. The difference between each 
spectrum and the initial spectrum at T = 10°C is shown 
in Figure [3}). The isosbestic point between these two re- 
gions is an indication that the decreasing peak reflects a 
decrease in the population of one species while the grow- 
ing peak is a sign of its conversion to a new species. This 
process is found to be associated with a decrease in the 
ASE peak, an indication that the population of DO 11 
molecules that generate ASE are being converted to ones 
that do not contribute to ASE. 

The positions of the two peaks in Figure [3)d are at dif- 
ferent energies than the corresponding peaks found dur- 
ing photodegradation, as shown in Figure [T];, suggesting 
that the product formed at higher temperature is not 
necessarily the damaged state. However, the most con- 
vincing evidence that an increase in temperature does not 
produce a damaged state is that the process is instanta- 
neously reversible, that is, the spectrum change follows 
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FIG. 3: (a) Absorbance spectrum of DOll molecules when 
doped in PMMA at various temperatures, (b) Difference in 
absorbance with respect to absorbance at 10°C. As the sam- 
ple temperature increases, the population of the bystander 
state increases as can be seen from the growing peak at 3.0eV 
while the population under main peak decreases. 
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FIG. 4: Height of the absorbance peak is plotted as a function 
of temperature. The data are fit to Equation [5] to get an 
estimate of the energy difference between molecular ground 
state and the bystander state, and it is found to be given by 
Eb = 73.4(±l.l)meF. 



the temperature with negligible delay. In contrast, re- 
covery of the damaged species takes several hours, even 
at elevated temperatures. We take this large difference 
in time scales as strong evidence that the bystander state 
is not a damaged species. 

The peak at higher temperature is not related to the 
decay product, and is most likely related to a different 
form of the DOll molecule, such as the DOll molecule 
in an excited vibronic state in the electronic ground state 
manifold, an isomer, a tautomer, or twisted intramolecu- 
lar charge transfer (TICT) state. The important point 
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is that whatever the species, it is not a decay product but 
does not contribute to ASE. Given that this product does 
not participate in the photodegradation process, we refer 
to it as a bystander state of the undamaged molecule. If 
the energy difference between the undamaged molecule 
and the bystander state is et, the temperature depen- 
dence of the undamaged population is given by, 

N 

^ l + exp[-..AT] 

Figure m shows a fit to the temperature dependence of the 
peak of the absorption spectrum shown in Figure [3] from 
which we get Sb = 73A{±l.l)meV. Note that at room 
temperature, the bystander population is about 5%. 

Generalizing Equation|3]to include the bystander state 
is somewhat complicated by the details of the damage 
mechanism. For example, does the pump beam damage 
both the DO 11 molecule in its undamaged and bystander 
state? Is the damaged state also characterized by two 
species? If the damage process removes molecules from 
the undamaged population, thermalization will result in 
an increase in the undamaged population from the reser- 
voir of molecules in the bystander state. Similarly, if a 
bystander state is associated with the damaged species, 
thermalization will result in a combination of both. Com- 
binations and permutations of these processes can lead 
to even more complex behavior. 

Because the population of the bystander state is small 
relative to the undamaged DOll population, we will 
generalize Equation 2] in the most straightforward way 
by simply replacing N with N\ thus excluding the by- 
stander state from consideration, or 



n[t;N'{T)J] 



{N' - al/l3) no 



no + {N' -no-iaI/l3))exp[-{N'l3-aI)t]' ^^'^^ 

By replacing N with N'^ we are not reducing the do- 
main size but the effective population that participates 
in decay and recovery of ASE. uq continues to refer to 
the initial undamaged population, but excludes DOll 
molecules in the bystander state. The theory can be 
easily generalized to account for complexities associated 
with individual systems. 

The theory embodied in Equation [10] is general in the 
sense that it applies to a broader set of systems than 
those described when motivating the derivation. For ex- 
ample, it has been suggested that the decay product 
might be a molecule with an expelled electron that is 
trapped in the polymer matrix. [15| an aggregate formed 
between two molecules, or between molecules that are 
bound upon the exchange of a proton, [llj All these mech- 
anisms can be equally represented by our model. 



D. Distribution of Domains 

The recovery process is based on groupings of 
molecules into generic domains and the dynamics were 



shown to depend on domain size. The domain size will 
be governed by the competition between attractive forces 
and thermal disordering, leading to a distribution of do- 
main sizes. Examples of systems with a distribution of 
domain sizes include ferromagnets . 1171 micellized surfac- 
tant solutions, [l^ liquid crvstals.(l9f and of course dye 
solutions. In each case, the equilibrium domain size dis- 
tribution is associated with a minimum of the free en- 
ergy. The domain size distribution can be derived in sev- 
eral wavs. j20l-l22j The most common and simple method 
is minimizing the Helmholtz free energy using a grand 
canonical partition function. We use this approach be- 
cause our system shares with these others the interaction 
between entities (in our case molecules, mediated by the 
polymer) that form a domain and thermal interaction 
that limit domain size. 

The partition function, z^, for a domain with N 
molecules is given by. 



Zff — exp 



A(A^- 1) 
kf 



= exp [7(iV - 1)] (11) 



where A is the free energy of a single molecule outside of a 
domain relative to a molecule within a domain such that 
the energy associated with a domain with N molecules 
is. En = —X{N — 1), k is the Boltzmann constant, T 
is the absolute temperature, and 7 = X/kT. When A is 
positive, energy is released when the molecule is added 
to a domain. The global partition function of a collection 
of domains is then given by. 



N 



(12) 



where f2Ar is the number of domains with N molecules. 
If our system has unit volume, JIat is the number density 
of domains of size N. 

The Helmholtz free energy, F, can be obtained from 
the partition function and simplified using Stirling's ap- 
proximation. 



F = -kTlnZ 



-kT 



N 



kTj2^N 
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ZN 



(13) 



The chemical potential of a domain of size N is therefore, 



dF 

an 



kT In- 



N 



N 



ZN 



(14) 



At equilibrium, the chemical potential of each 
molecule, whether it is a single molecule or part of a 
domain of any size, is the same. Equating hn /N and /xi, 
one obtains a relationship between the number density 
of domains of size N and the number density of single 
molecules as. 
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FIG. 5: Number density of domains Q{N; p, T) as a function 
of domain size at several temperatures (a) and concentra- 
tions (b). 



FIG. 6: Evolution of rj(N) as a function of time during con- 
tinuous pumping. 



by, 



The total number of molecules in the system is given 



N = l 



N 



J2 Ne-'^^-'^n^ 



N=l 



(16) 



where p can be interpreted as the average number density 
of molecules for fixed total volume, which is also simply 
related to the concentration of dye in the system. 
After some rearrangement of Equation 1161 we get 



(1 + 2^6-^)- Vl-H4peT 



(17) 



Substituting Equation 117! into Eguation ilSl we can write 
the number density of domains of size N, n{N) as a 
function of number density, p, 

m , o ^ / -I I ^ — -]N 



(1 + 2pz) - VI + 4pz' 



2pz 



N 



(18) 



where z = exp (■^). 

In the model, the density, p, and temperature, T, are 
experimentally controllable, with the free energy of a sin- 
gle molecule outside of a domain relative to a molecule 
within a domain. A, the only free parameter. Figure [5] 
shows a plot of the simulated distribution of U{N; p, T) at 
several concentrations and temperatures with the other 
parameters fixed. At higher temperature, the average 
domain size decreases as the thermal energy breaks the 
domains apart. 



E. The Full Integrated Model 

The theory of correlated molecules that describes the 
recovery process of a domain is governed by Equation 110! 
while the thermodynamic model describes the distribu- 
tion of domain size as governed by Equation [181 The 
effects of correlations and the statistical model for do- 
main size can be combined through an ensemble average 
to predict the number density of the undamaged popula- 
tion as a function of time, temperature, light exposure, 
initial non-equilibrium undamaged population, and con- 
centration of chromophores. 

In defining the domain size in the thermodynamic 
model, N refers to the total number of molecules in 
a domain, including the DOll molecule, the bystander 
species and the degraded species. However, the recovery 
process is hypothesized to be governed by interactions 
between DOll molecules with the bystander state re- 
moved, so the bystander species must be excluded. The 
mean number of undamaged molecules in a domain, n, is 
thus given by the ensemble average, 

n{t- p, T, I, no) = J2 WN; P, T) (19) 



n{t)n{N)dN (20) 
r]{N,I,t)dN, (21) 



where Equations [10] and [TH are used to evaluate the in- 
tegral in Equation!^ and N' is given by Equation [3| 

The distribution function r]{N, I, t) represents the pop- 
ulation of the undamaged species and is thus propor- 
tional to the contribution to optical absorbance of the 
undamaged species from each group of domains of size 
N . The total ASE emission from domain size N is also 
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FIG. 7: Evolution of ri{N) as a function of time during recov- 
ery after pumping for 10 minutes. 



FIG. 8: Simulated undamaged population decay as a function 
of time at several temperatures at particular dye concentra- 
tion and pump intensity. 



related to ri{N,I,t). ri{N)AN is the fraction of the un- 
damaged population living in domains of size between N 
and TV + AN. 

Figure[6]shows simulations of the evolution of ri{N, I,t) 
over time as the material is pumped. Area under 
each curve represents the ensemble average n(t). Since 
molecules in smaller domains degrade at a higher rate 
than molecules in larger domains, the peak in the dis- 
tribution shifts to the right upon photodegradation, as 
shown by the dashed curve. At long times, the distribu- 
tion function converges to an equilibrium shape as the 
decay and recovery rates balance each other in each do- 
main. When the pump is turned off, the damaged species 
in the larger domains will recover more quickly. As the 
larger domains recover, the smaller domains follow, re- 
sulting in the peak position shifting back to the left, as 
shown in Figure[71 The peak position as a function of do- 
main size during a run of decay and recovery thus traces 
a hysteresis loop. 

From Equation I21[ we can predict the behavior of the 
decay of dye-doped polymer as a function of dye con- 
centration and temperature. Figure [5] shows a series of 
simulated curves of the time dependence of ASE as pre- 
dicted for various temperatures. 



IV. TESTING THE THEORY 

Experimental details are described elsewhere [l^. It is 
important to note that due to spatial variations from spot 
to spot, it is difficult to get reproducible data even from 
a single sample. Typical point to point sample variations 
can yield as much as 30% variation in the signal, with 
larger variations from sample to sample. As such, most 
measurements require multiple runs to be averaged to 
decrease the effects of such variability. 

In determining the change of optical absorbance, the 



white light source - which focuses to a circular area, must 
be overlapped with the region that generates ASE. Since 
the pump light must be focused to a thin line to opti- 
mize ASE signal, the pump light and the wider white 
probe beam can never fully overlap. It is also difficult to 
assure that the overlap region is optimized, so, most of 
the probe might be missing the damaged area. To make 
matters worse, the pump beam can have hot spots that 
evolve over long data runs - which can take up to sev- 
eral days. Thus, some areas along the pump line may 
be fully damaged while others are close to pristine. This 
makes it difficult to definitively determine the change in 
the absorption spectrum for a damaged area. To over- 
come these issues requires painstaking adjustments and 
multiple runs in different spots on one sample as well as 
runs on different samples prepared in the same way. 
Empirically we find that the ASE intensity satisfies. 



Iase — 



(c/co)^ 



1 + ih/Ipump)P 



(22) 



and that the undamaged population is governed by the 
equation 



(23) 



for a particular pump inteirsity, where where Iq , p, cq and 
q are constants and where c is the dopant concentration 
and Ipump is the intensity of the pump. Figure [3] shows 
the ASE intensity as a function of the pump intensity for 
several concentrations of DOll dye as well as the fit to 
Equation 1221 The curves are normalized for the figure to 
make them overlap to show that they all have the same 
shape. The shape of the curve for the sample of 5g/l 
concentration is somewhat anomalous. 

Figure [TU] shows a plot of the ASE intensity as a 
function of dye concentration and a fit to the function 
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FIG. 9: ASE intensity plotted as a function of the pump in- 
tensity for samples of several concentrations. The inset shows 
the fit parameters 7o and exponent p determined from the 
data at each concentration. 
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FIG. 11: Normalized undamaged population (n) obtained 
from the ASE intensity is fitted to the correlation model to 
get the values of parameter p. 




TABLE I: Parameters determined for DOll in PMMA using 
an average pump intensity of Ip = 0.202W/cm^. 
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FIG. 10: Averaged ASE intensity at several dye concentra- 
tions. The data is fit to the Equation Iase oc (c/co)'' with 
Co = 5.0(±0.6) g/l and g = 2.6 ± 0.5. 



Iase oc (c/cq)'^, which yields cq = 5.0(±0.6) g/l and 
q = 2.6 ± 0.5. Thus, the ASE intensity can be used 
as a measure of the concentration of undamaged species 
through this relationship. 

As a test of the theory's ability to predict the depen- 
dence of the decay of ASE intensity as a function of time, 
samples of several concentrations are characterized using 
ASE at fixed temperature and pump intensity. The ASE 
intensity is converted to population of undamaged DOll 
molecules using Equation 1231 Initially the data from one 
concentration is fit to the model given by Equation [21] by 
independently varying all the parameters (a, /3, p, and A). 
Subsequently, the rest of the concentrations are fit to the 
model keeping constant the values obtained for a, /?, and 



Decay Rate 


I3(l0-*min-^) 
Recovery Rate Free 


\{eV) 
energy /molecule 


7.09(±0.13) 


3.22(±0.26) 


0.29(±0.01) 



A above and allowing only p, the number density of DOll 
molecules, to vary. A value of p is determined for each of 
the concentrations. 

Figure [Til shows representative data and fits to the the- 
ory for three concentrations. Six different concentrations 
are tested, and multiple runs at each concentration yield 
multiple values of the fit parameter p, which are aver- 
aged. The experimental uncertainty is determined from 
the spread in the data. Table [J summarizes the values of 
the three parameters obtained from the fits. These same 
values of the three parameters are applied to the the full 
set of data in this paper, a, j3 and A are properties of 
the material, and therefore should be constant for a given 
polymer and dye combination, though variations in the 
polymer due to materials processing may result in slight 
differences in these constants. 



Using the parameters determined from the decay data, 
the recovery data is predicted without the use of ad- 
justable parameters. Figure \\% shows the predicted pop- 
ulation as a function of time (curves) during recovery and 
the measured population of undamaged DOll molecules 
as determined from the ASE intensity. We note that dif- 
ferent samples were used to study recovery. Even so, the 
theory's prediction agrees with the data within experien- 
tial uncertainties for a range of concentrations. 

According to the model, p is proportional to the con- 
centration of dyes in the polymer. As such, we expect a 
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FIG. 12: Recovery of undamaged population (n) and a fit 
to the correlated chromophore model using the parameters 
obtained from fits to the decay of population during pumping. 



linear relation between p and concentration. Figure [T3l 
shows a linear fit to a plot of p as a function of the con- 
centration of a series of samples as determined from the 
amount of dye added to the polymer during its prepara- 
tion. The linear fit is consistent with the data. 



V. CONCLUSIONS 

We have presented a simple model for photodegrada- 
tion and recovery that depends on three parameters: the 
intensity-dependent decay rate, the recovery rate, and 
the free energy of a single molecule outside of a domain 
relative to a molecule within a domain. This model ac- 
counts for all observations of ASE and absorption spec- 
troscopy of DOll dye in PMMA as a function of time, 
pump intensity and concentration during decay and re- 
covery with one set of these three parameters. Further- 
more, the theory predicts the behavior as a function of 
temperature, which we are in the process of testing with 
experiment in which preliminary data is consistent with 
predictions. A more exhaustive experimental study will 
be presented in the future. 

We find that three parameters fully characterize a com- 
posite material made from a particular chromophore and 
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FIG. 13: Fit parameter p from decay fits as a function of dye 
concentration as determined during sample preparation. 



host polymer. As such, a, /3 and A - which should be 
calculable from first principles based on the underlying 
mechanisms - hold the key in providing a connection be- 
tween experiment and more fundamental theoretical con- 
siderations, a is related to the damage cross-section of 
the molecule, which is empirically found to decrease with 
increased domain size. /3, on the other hand, character- 
izes the self-healing process, with n/3 a measure of the 
collective strength of undamaged molecules to heal. The 
parameter A, on the other hand, determines the distri- 
bution of domain sizes, which are affected by tempera- 
ture and the concentration of the molecules in the sample 
when it is prepared. 

The theory is easily generalizable to account for other 
observations. For example, all studies to date where self 
healing is observed conclude that the degradation rate is 
proportional to the intensity, implying a linear mecha- 
nism of damage is responsible. For a nonlinear damage 
process. Equation [1] can be generalized by adding a non- 
linear function of the intensity. If an irreversible damage 
mechanisms acts along with a reversible one, an addi- 
tional species can be added to the model. 

The mechanism responsible for self-healing appears to 
be in the interactions between molecules that are some- 
how mediated by the polymer. Indirect evidence suggests 
that damage is associated with charge ejection from pris- 
tine molecules and the creation of a trapped charge [23j 
density in the polymer. [T5| We propose the hypothe- 
sis that in liquid solution, molecules that are damaged 
may break apart into fragments (perhaps charged) that 
by virtue of mixing in the liquid state takes them too 
far apart, on average, to recombine over a reasonable 
time frame. Once the molecular ions are neutralized by 
charges in the liquid through collisions and mass trans- 
port, reconstitution of the original molecule from the 
fragments is energetically unfavorable, and the process 
is thus irreversible. In this hypothesized mechanism. 



10 



the ionic fragments separate, but because of the poly- 
mer, remain in closer proximity to each other. Electro- 
static attraction between the fragments drives recombi- 
nation. If molecules are associated with others in a do- 
main, the larger fragments remain relatively stationary 
compared with the lighter ones. Thus, as higher num- 
bers of fragments are produced, attractive forces to the 
domain increases, enhancing self-healing, as experimen- 
tally observed. Dielectric screening due to the polymer 
may act to lengthen the healing time. 

Another possibility, though less likely, is that self heal- 
ing is a process that is a totally new phenomena analo- 
gous to Bose-Einstein condensation, which favors relax- 
ation into a Fock State in which all particles are in the 
same single-particle state. In this picture, the molecules 
are correlated due to positive exchange statistics of some 
sort, where a domain of undamaged molecules induces 
healing in proportion to their population. 

The source of correlations between molecules is not 
clear; but, there are many possibilities. It has been 
suggested that the association may be in the form of 
aggregation. jTli [23 | which can originate in either electro- 
static interactions, hydrogen bonding, or the formation 
of nanocrystalites. Alternatively, phonons in the polymer 
chains - which behave as bosons, may result in correla- 
tions. Or, the mechanism may be an altogether new phe- 
nomena. Any process in which an aggregate is of lower 
energy, or is more probable due to exchange statistics, is a 
viable candidate. Independent measurements are needed 
to sort them out. 

There are many experiments that can be brought to 
bare on the problem of mechanisms. The role of the 
polymer and its interaction with a guest molecule can 
be tested by varying the properties of each. For exam- 
ple, the fact that self healing is not observed in liquid 
monomer but in polymer suggests that there exists a crit- 



ical level of polymerization that leads to self healing. As 
such, one could measure the healing rate as a function 
of polymerization, is situ - while a dye-doped monomer 
is in the process of polymerizing - and correlate molec- 
ular weight with healing to determine if each domain is 
composed of molecules that are associated with a single 
polymer chain. Alternatively, the material can be heated 
through the polymer glass transition to determine if in- 
creased mobility of the polymer chains interferes with the 
healing phenomena. In addition, the polymer host and 
guest molecule can be changed to determine the impor- 
tance of structural and chemical properties to healing. 

Optical characterization, including various types of 
imaging, ASE, absorption spectroscopy, and fluorescence 
can be applied in tandem with other experiments such 
as photoconductivity to test a given hypothesis. For 
example, an observation of photoconductivity in coin- 
cidence with changes in population of damaged species 
as probed optically would support the hypothesis that 
charge ion generation plays a role. Suppression of heal- 
ing in the presence of a strong electric field would support 
the mechanism of charged molecular fragments, as would 
thermally stimulated discharge measurements. Finally, 
neutron scattering experiments could be used to directly 
probe correlations and aggregation. 

In summary, the model that we present here predicts 
all of the observations in terms of three parameters and 
is based on the hypothesis that self healing originates 
from a collective phenomena in which pristine molecules 
induce healing in a damaged molecule in proportion to 
domain size. Similarly, the decay rate decreases with 
domain size. Future planned experiments will be used to 
zero in on the mechanisms. 
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